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Taking the sum rule for the p meson as an example, we study the possibility 



^ ■ of explaining the phenomenological enhancement factors for certain terms in the 



vacuum expectation value of the operator product expansion in the QCD sum rule. 
We take a QCD motivated extended Nambu-Jona-Lasinio (ENJL) model as the low 
energy effective Lagrangian for QCD with which we calculate the nonperpturbative 
contributions to the vaccum condensate expansion to obtain the enhancement factors. 
Our result shows that such nonperturbative contrbutions can cause large enough 
enhancement factors which can be consistent with the phenomenological values. 
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I. INTRODUCTION 

QCD sum rule [I] is a useful semi-phenomenological approach to hadron physics in 
which nonperturbative vacuum condensate (VC) effects are taken into account. The sum 
rule relates physical observables to certain VC's through the vacuum expectation values 
(VEV's) of products of current operators (Tj A (x)j B (0)) with the vacuum condensate 
expansion (VCE) 

• / d 4 x e** (Tj A (x)j B (0)) = F AB (q) £ C n (Q 2 ) (O n ) , (1) 

where F AB (q) is a tensor factor characterizing the tensor structure in the indices A and B , 
Q 2 = —q 2 , (O n ) 's are VC's of various dimensions. Usually, up to the VC's of dimension-6 
operators (four-quark operators) are taken into account in (1). In most applications of the 
QCD sum rule, the following approaches are taken: 

(i) as a basic assumption, the coefficients C n 's, at least in the first few terms, are assumed to 
be just the Wilson coefficients in the standard operator product expansion (OPE) determined 
merely from perturbative QCD (PQCD), and nonperturbative effects all reside in the VC's 
(On) 's 0; 

(ii) (O n ) 's are treated as free parameters determined by experimental inputs, and to reduce 
the number of free parameters , people often make the simple factorization approximation 
to express the four-quark condensate (4QC) in terms of the two-quark condensate (2QC) in 
the following way []]]. 

^r^VT^) « N- 2 [ TrT 1 Trr 2 - Tr^a) ](^) 2 , (2) 

where ip is the three-flavor light quark field, the Lorentz structure of r 12 is 1 for scalar 
(S), for pseudoscalar (P), 7 M for vector (V) , 757^ for axial- vector (A) , a av for tensor 
(T) (for color and flavor non-singlet, the color and flavor group generators A a /2 , should 
also be included ), and iV = 36 is a normalization constant JI||. 
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However, a careful analysis of the QCD sum rule for the p meson shows that 
the theoretical results are not consistent with the experiments and phenomenologically the 
demension-6 4QC term should be enhanced by a factor 

k 4 q ~ 3 - 6.6 (3) 

for fitting the data ||. Studies of other processes with the same approach such as the 
sum rules for baryons fU and pseudoscalar meson 0] lead to similar conclusion that an 
enhancement factor around k^q « 4 is needed for fitting the data. 

There is an alternative analysis from the finite energy sum rule which leads to similar but 
somewhat different results that the dimension-6 4QC term and the dimension-4 gluon con- 
densate plus two-quark condensate terms (Actually, the gluon condensate term dominates.) 
should be enhanced by factors 

« 5 - 8, kq ~ 2 5 , (4) 

respectively [J5J. 

So far there is no successful theoretical explanation of the enhancement factors. In our 
previous paper ||, we kept assumption (i) and looked at the possibility of explaining the 
enhancement factor k^q in (3) by taking account of the non-factorized parts of the 4QC's, 
i.e. modifying approximation (ii). We took the QCD motivated Nambu- Jona-Lasinio (NJL) 
model by Bijnens, Bruno, and de Rafael |7j as the low energy effective Lagrangian of QCD to 
calculate the non-factorized parts of the 4QC's by means of a nonperturbative method for the 
effective potential for local composite operators ||, and the obtained non-factorized parts 
of the 4QC's are of the next-to-the-leading order in the 1/N C expansion. Our conclusion 
is that for reasonable values of the parameters in the model, the obtained enhancement 
factor k 4 q is only a few percent larger than unity, so that it is too small to account for 
(3). The next reasonable trial is to consider certain modification of assumption (i) since the 
p meson lies below the chiral symmetry breaking scale wherein nonperturbative effects are 
important. Although an argument considering nonperturbative effects from instantons in 
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the dilute gas approximation has been given in Ref. [0 expecting that assumption (i) may 
work up to about ten terms in (1), the analyses in Refs. @ |J show that the real situation 
is not so optimistic. In this paper, inspired by a recent paper by Yamawaki and Zakharov 
0, we study the possibility of explaining the enhancement factors by taking into account 
certain nonperturbative contributions to the VCE in which nonperturbative effects not only 
reside in the condensates but also affect the coefficients in the VCE so that assumption 
(i) is modified. We take the whole extended Nambu- Jona-Lasinio (ENJL) Lagrangian in 
Ref. as the low energy effective Lagrangian of the quark and gluon system, and calculate 
the nonperturbative contributions to the VCE from the NJL four-fermion interactions to 
the precision of the leading order of 1/N C expansion. Our results show that, taking 
only four-fermion interactions of the scalar and pseudoscalar types ]7), such modification 
mainly affects the 4QC terms in the VCE by causing additional 4QC's with different tensor 
structures which leads to a fairly large contribution to k^q . For reasonable parameters in 
the model, it is really possible to obtain a large enough k±q consistent with (3). We shall 
also show that, when including certain vector and axial- vector four-fermion interactions in 
the ENJL Lagrangian |7|], it is also possible to explain the enhancement factors (4). This 
means that for low energy processes, QCD sum rule may give better results if the original 
assumption (i) is modified by taking account of more nonperturbative QCD effects. 

This paper is organized as follows. Sec. II is the theoretical aspect of this study. In Sec. Ill 
we present the calculations of the nonperturbative contributions to the VCE up to the 4QC 
terms from the ENJL model and the explanation of the enhancement factors. Conclusions 
are given in Sec. IV. 

II. THE THEORETICAL STRATEGY 

To study nonperturbative contributions, it is convenient to take a low energy effective 
Lagrangian of QCD including nonperturbative effects. So far there is no successful low energy 
effective Lagrangian derived from the first principles of QCD. The QCD motivated ENJL 
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model by Bijnens, Bruno and de Rafael contains reasonable QCD ingredient including 
chiral symmetry breaking, and can lead to rather successful phenomenological results [0. 
So we take this model as the low energy effective Lagrangian for our calculations. The 
Lagrangian in the ENJL model in Ref. |7j is 

Cqcd = Cqcd + £ S njl + £njl > ( 5 ) 

where A x is a momentum cut-off below which the model serves as the low energy effective 
Lagrangian for QCD, Cq ) cd is of the same form of the QCD Lagrangian for momentum 



below A x , and 



iVcA X ab 

[(^) 2 + (V^) 2 ] (6) 



3iV c A x 
Air 2 C 8 

iV c A x i= i 



SN C A 2 X 
8n 2 G v 



[(^VO 2 + (^7s7^) 2 ] (7) 

v x3 E[(^^) 2 + {^l^f , 
in which and Gy are two coupling constants treated as free parameters, and the flavor 
group generator £ 4 is normalized as tr(tjt,) = |5y . It is argued in Ref. that the NJL 
type Lagrangian C S pf JL + C V N ^ L may be understood as coming from integrating out the 
high momentum modes of quarks and gluons in the fundamental theory of QCD (However, 
Gs and Gy are taken as two independent free parameters here.) and is regarded as the 
part of Cqcd responsible for chiral symmetry breaking, and Cqcd provides perturbative 
corrections to the broken chiral symmetry state. This supports the idea of the chiral quark 



model |0|. 

Our strategy is that we apply the Lagrangian (5)-(7) to energy scale lower than A x , and 
take the ordinary PQCD approach for the energy range higher than A x where nonperturba- 
tive effects are not important. Since the structure of Cq^d is the same as the original QCD 
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Lagrangian, the pure perturbation results in the present appraoch are the same as those 
from the ordinary PQCD. The NJL Lagrangian Cff JL + C^ A JL not only gives rise to chiral 
symmetry breaking, but also give nontrivial contributions to the VCE which is what we are 
going to calculate. 

There are three independent free parameters G$ , Gy , and A x in the ENJL model 
in Ref. 0. They are related to the dynamical quark mass Mq , the quark axial- vector 
coupling constant qa , and the vector meson mass M p by the following relations f?J 

1/G S = (M Q /A x ) 2 r(-1, (M Q /A X ) 2 )(1 + 7 _0 , (8) 

9A = 1 + 4G V (M Q /A X )*T(0, (M Q /A X )*)(1 + 701 ) ' (9) 

A 2 = ^M p 2 G y r(0, (M Q /A X ) 2 )(1 + 703 ) , (10) 

where T(n — 2, (Mq/A x ) 2 ) is the incomplete gamma function 0, and 7_i , 701 , and 
703 are perturbative corrections. Actually Mq , A x , and are taken as the input 
parameters to fit the data in Ref. 0. In this paper, we calculate the contributions of the 
NJL La grangian jCfJj L + to the precision of the leading order in the l/A^-expansion, 
i.e. the uncertainty of the calculation is of the order of 30%. With this uncertianty, we shall 
neglect the small quatities 7 _i, 701 and 703 throughout this paper. Moreover, we shall also 
neglect the non-factorized parts of the 4QC's given in Ref. M since they are of higher order 
in the expansion. 

In calculating Feynman diagrams to the precision of the leading term in l/iV c -expansion, 
we should sum up all open-end chain diagrams attached to each quark line (cf. Fig.l). 
Note that only the flavor-singlet scalar four- fermion interaction contributes to these chain 
diagrams. It is easy to see that the effect of summing up such chain diagrams attached to 
a quark line is just providing a dynamical mass term Mq (arising from (tptp) 7^ ) added 
to the current quark mass m to form the full quark mass of this quark line. The relation 

r> 
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Fig. 1. Chain diagrams in the calculation of the leading contribution in 1/N C expansion. 



between Mq and (ipip) can be obtained as follows. In the large- N c limit, the gap equation 

m = -^Kn-i,(M 2 Q /A 2 x )) , (ii) 

in which the chiral limit relation (uu) = (dd) = (ss) = ^(4>ip) is used, the function 
r(— 1, (Mq/A^)) can be eliminated by using (8). neglecting 7_ 1 in (8), we get from (11) 
and (8) 

M " = -^f<« • (12) 

Hence Mq is proportional to (ipip) ■ Since QCD sum rule concerns the VCE in which 
various dimensional quark condensates are put as explicit terms, all the effects of Mq 
should be put into the explicit quark condensates in the VCE by definition. Considering 
these nonperturbative effects, (1) becomes 

i I d*xe**{j A (x)j B {0)) = F AB (q)Y,C n (Q 2 )(O n ) , (13) 

J n 

in which the new coefficient C n (Q 2 ) in the resulting VCE is not the same as the original 
Wilson coefficient C n (Q 2 ) in the standard OPE. C n (Q 2 ) contains nonperturbative effects 
from the NJL Lagrangian. Generally, C n (Q 2 ) is a function of A x , Gs, Gy, and the 
current quark mass m . In the theory with the Lagrangian (5)-(7) ; it is (13) rather than (1) 
that really matters in the QCD sum rule. Note that when calculating C n (Q 2 ) , the quark 
mass in every quark line in the Feynman diagrams should be understood as the current quark 
mass m . Since the current quark mass m is much smaller than Q 2 ~ M 2 , we treat 
m 2 /Q 2 as a small perturbation in the following calculations. 
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III. CALCULATION OF C n 's AND THE EXPLANATION OF k 's 

In the sum rule for the p meson, the currents j A , j B are taken to be vector currents 

U = ^(«7m w ~ hud) = 4>t 3 ?p , (14) 

where £3 is the flavor group generator in (6)-(7). The VCE of the product of the currents 
takes the form 

i J dx (0|Tj M (xH,(0)|0> = {q,q v - q 2 g, u )Tl(q 2 ) , 

n(Q 2 ) = Cj(i} + c 2Q $M u , d ip) + c G (G%G a n + c 4Q 53(^r^r^> + • • • , (is) 

r 

where M u ^ is the u—, d— quark mass matrix. We are going to calculate the four terms 
on the right-hand-side of (15) including the nonperturbative contributions from the NJL 
Lagrangian in (5)-(7). 



1. Model with G s ^ 0, G v = 0. 

There are five different ways of fitting the data presented in Ref. 0, which determine 
different sets of values of Mq , A x , and qa , and the predictions are all successful. 
The simplest one of the fits is their Fit 4 in which = 1 ■ From (9) we see that this 
corresponds to Gy = which makes the calculation of the nonperturbative contributions 



easiest. Moreover, it is explained by Weinberg JL1] that, in the large iV c limit, should 
actually be unity in the constituent quark model. In view of this, we take the set Fit 4 in Ref. 
in this subsection, and by means of (8)-(10), it leads to Gs = 1.19, Gy = 0, A x = 667 
MeV. Since the uncertainty in our calculation is of the order of 30%, we are not taking these 



in 




1 1 



Fig. 2. Example of Feynman diagrams for leading contributions to 
Ci(I) . (a) PQCD contribution; (b) Cf/jL contribution. 



precise values seriously and, instead, we shall consider a reasonable range of the parameters. 

Now we calculate the coefficients in (15). Examples of the Feynman diagram of the 
leading contribution of PQCD and jCff JL to Cj(I) is shown in Fig. 2. This has been 
studied in Ref. ||. The first nonvanishing contribution comes from the term suppressed by 
Gsm 2 /Q 2 which is negligibly small compared with our theoretical uncertainty. Thus In the 
present appraoch Cj(Q 2 ) is fairly well given by the standard Wilson coefficient Cj(Q 2 ) 
from PQCD given in Ref. M, i.e. 



The PQCD result of the C 2 q{iPMi(j) term is proportional to the current quark mass 
m |J [cf. Fig. 3 (a)]. Examples of the Feynman diagram for the Cjf JL contribution to the 
C^Q^Mij)) term is shown in Fig. 3(b) (This corresponds to the effect linear in Mq .). It is 



negligible relative to the PQCD contribution to the present precision. Hence C2q{Q 2 ) i n 
the present approach is also mainly given by the standard Wilson coefficient C2q{Q 2 ) from 




i , Q 2 



(16) 



easy to see that the first nonvanishing contribution is of the order of m— — and is thus also 

1)2 



PQCD given in Ref. IE, i.e. 



1 



(m u (uu) + m d {dd)) . 



(17) 
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Fig. 3. Example of Feynman diagrams for leading contributions to C2Q{ipM u dip) ■ (a) PQCD 

contribution; (b) Cff JL contribution. 

The situation of the Ca{G a G atlv ) is similar. Examples of the Feynman diagrams for 
leading PQCD and jCff^ L contributions to this term are shown in Fig.4. The first nonva- 
nishing contribution in Fig. 4(b) is also of order Gsm 2 /Q 2 which highly suppresses the 
nonperturbative contribution. Thus Cg{Q 2 ) is also maily given by the standard Wilson 
coefficient Cg(Q 2 ) from PQCD given in Ref. PJ, i.e. 

~ rv 

C G (G« u G a n « C G {G%G*n = ^Q- 4 (G°„G a n . (18) 

For the 4QC term, the situation is completely different. Examples of the Feynman 
diagrams for leading order PQCD and CffjL contributions are shown in Fig.5. Unlike the 
previous cases (cf. Figs. 2(b), 3(b), 4(b), Fig.5(b) is of tree-level rather than loop-level, 
therefore this contribution is large. The leading PQCD [Fig.5(a)] formula given in Ref. 
can be written as 



r 

2ixa 



PQCD V 



Q 



6 



(19) 



where C\q(Q 2 ) is the standard Wilson coefficient. Note that there are only VEV's of 
products of vector and axial-vector currents. The calculation of Fig.5 (b) is straightforward 
but lengthy. There is also contribution to the 4QC term from the expansion of (ijj(x)ijj(O)) 
arising in Fig.3(b). Adding the two contributions together, we get 



IF) 



Ifi 
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Fig. 4. Example of Feynman diagrams for leading contributions to CciG^G 11 ^) . (a) PQCD 



contribution; (b) £jvjl contribution. 



2tt 2 G'c 



NJL 



N C A 2 X Q* 



(20) 



where t± = ti ± ^2 • We see that (20) contains the VEV's of products of scalar and 
pseudoscalar 4QC terms which are not suppressed by the current quark mass m , and thus 
are not small corrections to (19). Adding (19) and (20) together, we get the total 4QC term 



PQCD 



(21) 



NJL 



Compared with the standard OPE formula (19) used in Ref. 0, (21) contains extra significant 
terms which enhance the 4QC term in the VCE (15). 

The above results are consistent with the analysis in Ref. 0, i.e. the first three terms in 
(15) are as usual, while the 4QC term is enhanced. With (21), we can define the enhancement 
factor k 4 q in this approach as 



K 4Q 



[22) 



where the numerator is given by (21) and the denominator is given by (19). The 4QC's in 
(19) and (21) contain their factorized parts given by (2) and non-factorized parts given in 



IS 



19 



Fig .5 Example of Feynman diagrams for leading contributions to X^r (tpTipipTip) 



. (a) PQCD contribution; (b) Cff JL contribution. 



Ref. ||. As we have stated, to the present precision we neglect the non-factorized parts. 
The factorized part of (19) reads || 



r 



pqcd_ si I w 



(wr- (23) 



factorized 

The factorized part of (20) obtained from (2) is 

4tt 2 Gc 



c 



njl . 9N C AIQ 4 

j actorized A 



<^> 2 - (24) 



From (12) we see that these factorized parts correspond to Mq contributions to the 4QC 
term. The explicit formula for k^q for a given Q 2 in the present approach is 

187rAy; s Q 2 

»«» » 1 + 7 (Ag - 1KA| • < 25) 

We now estimate the value of k 4 q in (25). We take iV c = 3 and the PQCD coupling 
constant a s (M p ) « 0.48 which is obtained from the two-loop evolution formula from the 
world average value a s (Mz) = 0.118 [|l^]. From (8) we see that G,s and Mq/A 2 , are 
related. We take Mq = 199 MeV from Fit 4 in Ref. 0, and take A x around the chiral 
symmetry breaking scale A x ~ 4^/^ ~ 1.2 GeV. The values of k^q with Q 2 = M 2 for 
various values of A x are shown in Table 1. 
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Table 1 Values of k^q with Q 2 = M 2 for various given values of A x in the 

G s + 0, Gy = model. 



A x (GeV) 


G s 


K 4 Q 


1.25 


1.10 


3.7 


1.10 


1.15 


4.7 


0.90 


1.20 


6.7 


0.80 


1.25 


8.3 



We see that, with the uncertainty of the present calculation, the values of k 4 q are consistent 
with the values in (3). Actually in the analysis in Ref. 0, Borel transformation 

o, lim r^iv^ Q2 y i ^^) n ( 26 ) 

Q 2 /n^M 2 

was applied. We have also looked at the result after applying the Borel transformation (26) 
with M « M p . The results of k 4 q is larger than those listed in Table 1. For example, for 
A x = 1.25 GeV and 1.10 GeV, we have k^q = 6.5 and 8.3, respectively. These are still 
consistent with the phenomenological value in (3) considering the uncertainty of the present 
calculation. Thus the present theory with Gs ^ and Gy = is compatible with the 
analysis in Ref. || that the first three terms in the VCE (15) are the same as those based 
on the basic assumption (i), while the 4QC term should be enhanced by a factor k^q given 
in (3). 



2. Model with G s ^ 0, Gy ^ 0. 

The above model with Gs ^ 0, Gy = is not compatible with the enhancement 
factors (4) given in Ref. || in which k g is not unity. Since the £^-j L interactions give 
significant contributions to the first three terms in (15), it is interesting to see if a model 
with Gs 7^ and Gy ^ can be compatible with (4). 

We first consider the £]vjz, contributions to the 4QC term. The Feynman diagam for the 
leading contribution is the same as Fig. 5(b) with C^f JL replaced by C^j L . After straight- 

77 



forward calculations, we find that the factorized parts of ((^ ^t^ip) 2 ) and ((■07^75^3^') 2 ) 
cancel each other so that the net contribution of the C^j A JL interaction to the 4QC term is 
just the corrections from the non-factorized parts which is negligible. So that the value of 
k^q is not significantly affected by Gy 7^ . Let us take Fit 1 in Ref. J7J as an example in 
which g A = 0.61 , A x = 1.16 GeV, M Q = 265 MeV (7^ = 701 = 703 = in Fit 1). Fixing 
this Mq, the values of Gs and k^q with Q 2 = M 2 for various values of A x are listed in 
Table 2. We see that these values are consistent with (4) considering the uncertainty of the 
present approach. 

Table 2 Values of k^q and k g with Q 2 = M 2 for various given values of A x in the 

G s ^0, G v / model. 



A x (GeV) 


G s 


Gy 


K 4Q 




1.16 


1.22 


1.26 


4.4 


2.5 


1.00 


1.28 


1.06 


5.8 


2.7 


0.90 


1.40 


0.94 


7.5 


2.9 


0.80 


1.42 


0.84 


9.3 


3.1 



We then look at the £]vjl contributions to the gluon condensate term (the two-quark 
condensate term is relatively much smaller). The Feynman diagram is still Fig.4(b) with 
Cjfj L replaced by £]vjl- ft * s eas y to see that only the (ip r f^t 3 ip) 2 term in (7) contributes. 
This gives rise to 



C G (Q 2 ) = C G (Q 2 



+ C G {Q 2 



PQCD 

Thus we have the enhancement factor 



C G (Q 2 ) 
k g = ^ 



NJL 



1 + 



4GVQ 2 , Q 



N C A 2 X 



In 



C G (Q 2 



(27) 



C G (Q' 



l + 4^1nQ 2 



iV c A x 



In the QCD sum rule, the scale n is supposed to be the inverse of the confinement radius, 



R conf 0- TllUS P 



y MS 



280 MeV. The values of Gy and k g with Q 2 = M 2 for 
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various values of A x are listed in Table 2. We see that these are also consistent with the 
k g in (4). 

The £)vjl interactions will also give rise to the same enhancement factor as kq for 
the (7/(1) term For Q 2 « M 2 > /i 2 , the PQCD value of Cj [cf. (16)] is negative. 
Therefore this nonperturbative contribution makes Cj more negative. We expect a global 
analysis in the spirit of Ref. || including the Ci(I) term and test the present prediction. 

IV. CONCLUSIONS 

In this paper, taking the p-meson sum rule as an example, we have tried to modify the 
basic assumption of the QCD sum rule that the coefficients C n (Q 2 ) 's in (1) are merely 
determined by PQCD, and studied the possibility of explaining the phenomenological en- 
hancement factors for certain terms in (1) || || by taking account of nonperturbative 
contributions to the VCE which lead to (15). We take the QCD motivated ENJL model in 
Ref. |7[ as the low energy effective Lagrangian for QCD, and have calculated the nonpertur- 
bative contributions from the NJL Lagrangian to obtain the new coefficients C n 's to the 
precision of the leading order in the 1/N C expansion. 

It is interesting that taking Fit 4 of Ref. f7|, ( Gs ^ 0, Gy = 0), our result shows that the 
first three terms in (15) are not affected much by the nonperturbative interactions, while the 
4QC term contains extra terms with new tensor structures coming from the nonperturbative 
interactions and they are large compared with the PQCD result. Regarding these extra 
contributions as the theoretical source of the enhancement factor k 4 q , the obtained k^q 
[cf. Table 1) is consistent with the phenomenological value (3), considering the theoretical 
uncertainty in the present calculation, for reasonable range of the parameters in the ENJL 
model. 

Taking Fit 1 of Ref. J7|, i.e. Gs ^ 0, Gy ^ 0, we have got another result that the 
obtained k^q and k g listed in Table 2 are all consistent with the results (4) given in Ref. 
||. However, in this model, Ci(I) will also be larger than Ci(I) by the same factor 
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as listed in Table 2, but this has not been discussed in Ref. ||. We expect a further 
analysis of the Cj(I) term in the spirit of Ref. || to test whether the present prediction 
by Fit 1 is really right. 

We conclude that it is possible to explain the enhancement factors in both Ref. || and 
[|5| by taking different sets of parameters in the ENJL model in Ref. J7|. Within the 30% 
uncertainty of the present calculation, the model with Gs ^ 0, Gy = can give a 
reasonable explanation of the analysis in Ref. with the value of k^q given in (3), and the 
model with Gs ^ 0, Gy ^ ENJL model [0] can possibly explain the enhancement factors 
(4) given in Ref. || . In the present approach, the physical origin of the enhancement factors 
is essentially the nonperturbative contributions to the VCE which makes the VCE different 
from the original assumption (i), whereas the non- factorized part of the 4QC given in Ref. 
H is only a minor source of the enhancement factors. More precise nonperturbative QCD 
approach along this line is expected for developing practically useful algorithm to improve 
the application of the QCD sum rule to low energy processes. 
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